Monte Carlo generators for top quark physics at the LHC 



B.S. ACHARYA(l), F. CAVALLARI(2), G. C0RCELLA(3)(4)(*), ^ j^j SiPio(5) and 

G. Petrucciani("*) 

C') Abdul Salam International Center for Theoretical Physics and INFN, Sezione di Trieste, 
Italy 

(^) INFN, Sezione di Roma, Italy 

(^) Museo Storico della Fisica e Centra Studi e Ricerche E. Fermi, Roma, Italy 
(*) Scuola Normale Superiore and INFN, Sezione di Pisa, Italy 
(^) Universitd di Bologna and INFN, Sezione di Bologna, Italy 



Summary. — We review the main features of Monte Carlo generators for top 
quark phenomenology and present some results for ti and single-top signals and 
backgrounds at the LHC. 

PACS 14. 65. Ha - Top quarks. 

PACS 12.38.Bx - Perturbative calculations. 



1. Introduction 

In order to perform precise measurements of properties of top quarks at the LHC, 
the use of reliable tools will be essential. Extensive work has been carried out to im- 
prove parton-shower and matrix-element Monte Carlo (MC) codes for tt and single-top 
signals and backgrounds. Parton shower generators, such as HERWIG [1] or PYTHIA 
[2], simulate top production and decay using leading-order (LO) matrix elements and 
describe multiple radiation in the soft or coUinear approximation. HERWIG showers 
satisfy angular ordering, valid in the soft approximation after azimuthal averaging, while 
PYTHIA orders its cascades according to transverse momentum or virtuality, with an 
option to veto non-angular-ordered emissions. To simulate hard and large-angle radia- 
tion, both PYTHIA and HERWIG have been provided with matrix-element corrections. 
HERWIG splits the physical phase space into a soft/coUinear region, where one trusts 
parton showers, and a region, corresponding to hard and large-angle radiation, where 
the exact amplitude is used. Moreover, HERWIG corrects the 'hardest-so-far' emission 
in the shower. The PYTHIA implementation is somewhat different: the parton shower 
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approximation is used throughout the entire physical phase space and the first emission 
is generated exactly. HERWIG [3, 4] and PYTHIA [5] include matrix-element correc- 
tions to top-quark decay, but not to top production. Such corrections are also applied to 
W-boson production, which is matched to the -I- 1 jet process [6, 7], one of the back- 
grounds for top production. At the end of the parton cascade, HERWIG and PYTHIA 
simulate the hadronization transition according to the cluster [8] and string [9] models, 
respectively. As for the underlying event, both programs deal with it by allowing multi- 
ple interactions [10, 11]. For this purpose, HERWIG is interfaced with the JIMMY code. 
HERWIG++ [12] and PYTHIA 8 [13] are the new object-oriented versions of the two 
codes, written in C-l — h. PYTHIA 8 implements only transverse-momentum evolution, 
with initial/final-state radiation and multiple interactions interweaved in a common or- 
dering. HERWIG-I— I- includes a mass-dependent term in the splitting functions, leading 
to a better treatment of the radiation off heavy quarks and heavy-hadron fragmentation. 

Besides HERWIG and PYTHIA, other parton showers generators are available. The 
ARIADNE [14] code uses PYTHIA hard scattering and hadronization, but its showers 
are dipole cascades, i.e. 2 — > 3 branchings in the soft limit. The evolution variable 
of ARIADNE is transverse momentum: in fact, implementing dipole emissions, there 
is no need for azimuthal averaging and angular ordering. The SHERPA [15] generator 
implements showers which, like in PYTHIA, are ordered in virtuality, with an option to 
reject non-angular-ordered emissions. Hadronization occurs following the string model. 
A feature of ARIADNE and SHERPA is that matrix-element matching is implemented 
according to the CKKW-L [16, 17] method. Jets are clustered according to the kx- 
algorithm [18], with a threshold ycut on the fcr-clustering variable yij. Exact matrix 
elements, weighted by the Sudakov form factors, and parton-shower approximations are 
then used above and below ycut- 

The programs mentioned so far, even after matrix-element matching, still yield LO 
rates. Next-to-leading order (NLO) cross sections and observables are instead given 
by the codes MCQNLO [19] and POWHEG [20], both available for top production. 
MC@NLO implements the hard scattering at NLO and employs HERWIG for show- 
ers and hadronization: its drawbacks are its shower-model (HERWIG) dependence and 
its simulation of events with negative weights. Such problems have been overcome 
by POWHEG, which can be interfaced to any shower model and does not have any 
negatively-weighted event. The first POWHEG emission is the hardest in transverse mo- 
mentum {pt) and is generated exactly at NLO; the subsequent cascade is ordered in pT- 
To include angular ordering, one would need truncated showers, whose implementation 
is in progress. 

Finally, in order to simulate multi-jet final states in top-quark events and backgrounds, 
such as W/Z+ jets, the best available tools are the so-called 'matrix-element genera- 
tors', such as ALPGEN [21], MadGraph/MadEvent [22], CompHep [23], HELAC [24], 
GR@PPA [25] and WHIZARD [26], which include LO multi-parton amphtudes. Such 
programs can be interfaced to HERWIG or PYTHIA for showering and hadronization, 
according to the Les Houches accord. The ALPGEN code, for example, includes the 
following processes {V = W oi Z): tt + 6 jets, V + 6 jets, Wbb + 6 jets, VV + 3 jets. It 
uses the so-called MLM prescription to match its partons with the subsequent showers. 
Roughly speaking, one starts, e.g., with W + n partons and defines a cone jet algorithm. 
After showering, a fc-jet event is accepted only if k > n and each of the n original partons 
is clustered into a different jet. 
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Fig. 1.: Transverse momentum of the tt pair at the LHC, according to HERWIG, 
MC@NLO and the NLO calculation (left), and according to POWHEG and MC@NLO, 
for stable top quarks (right). 



, , , 


... 1 .... 1 .... 1 .... 1 . 

LHC 


0.25 




d(7/dy, (pb)- 








0.20 






" 0.15 




Y Solid: ALPGEN, Sj 


■S 0.10 

\ 




' M: MC^NI.O T^x - 


b 






0.05 
















0.00 



dashes: MC@NLO 
solid: POWHEG 




-4-2 2 4 

y of hardest jet 



Fig. 2.: Rapidity of the hardest jet in tt{g) events at the LHC according to ALPGEN and 
MC@NLO (left), and at the Tevatron according to POWHEG and MC@NLO (right). 



2. Results for top-quark signals and backgrounds 

We would like to present some results for top-quark signals and backgrounds at the 
LHC, using the programs described above. Following [27], in Fig. 1 we present the 
transverse momentum of the tt pair at the LHC according to HERWIG, MC@NLO and 
the parton-level NLO calculation, with all distributions normalized to the same area. At 
small p^*^, where one is mostly sensitive to soft/coUinear parton radiation, MC@NLO 
and HERWIG agree, while the pure NLO calculation is far above the two codes. At 
large p^*\ it is instead the NLO computation which is reliable: MC@NLO and the 
NLO prediction agree, while HERWIG underestimates the large-p^*^ rate. In Fig. 1 
we also compare MC@NLO and POWHEG for stable top-quark production: up to a 
small discrepancy at very low p^*^ , the agreement between the two codes is remarkable 
[28]. In Fig. 2 we present the rapidity of the hardest jet yi, yielded by MC@NLO and 
ALPGEN: we have clear disagreement around yi = 0, with ALPGEN giving more events 
than MC@NLO. As pointed out in [29], this is due to the fact that MCQNLO does not 
generate enough hard scatterings with ttg in the final state. This is confirmed by the fact 
that POWHEG, whose first emission is always generated at NLO, yields a j/i-spectrum 
similar to ALPGEN (see Fig. 2). 
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Fig. 3.: Left: relative transverse momentum of the hardest-jet partons with respect to 
the jet axis for single-top production at the Tevatron. Right: lepton pT for W -\-2 jets 
with a h tag, with the W decaying leptonically, at the LHC. The highest cross section is 
yielded by the subprocesses (solid), with V766[-hj] (solid) and W{hh)j (dashes) 

giving similar contributions. The NLO/LO ratio is higher for 
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Fig. 4.: Left: uncertainty Amt due to the underlying event in PYTHIA. The bands 
to the Amt > side account for jet-energy-scale corrections. In both cases, results on 
the left- (right-) hand side correspond to evolution in transverse momentum (virtuality). 
Right: Dependence of the tt cross section on mt for few values of renormalization and 
factorization scales. The lines on the top correspond to ~ '7it/2, in the middle to 
= mt, on the bottom to /.i/? = 2mt. 



For single-top production, a comparison among some of the available codes was per- 
formed in [30] . Fig. 3 shows the total transverse momentum of all partons in the hardest 
jet, relative to the jet axis, using HERWIG, MC@NLO and the NLO calculation. The 
agreement between HERWIG and MC@NLO is acceptable over the full range: in fact, 
p^j.gi is a sufficiently inclusive observable, so that NLO effects mainly result in an overall 
i^T-factor. The NLO result is sharply peaked at p^^ci = 0, since at NLO a jet often coin- 
cides with a single parton. At large p^rei^ the NLO result overestimates the prediction 
of the two MC codes, where a large-p^j.^,) hadron likely belongs to another jet. Among 
the backgrounds to single-top production, W + 2 jets, at least one with a b quark, were 
studied in [32], making use of MCFM [31], a NLO parton-level code. Fig. 3 shows the 
Pt distribution of the lepton from W decay, for the Wbb + X, Wbj + X and W{bb)j 
processes. At the LHC, Wbj + X is dominant, while the other two backgrounds are 
roughly comparable. 
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Fig. 5.: Left: S-liadron spectrum in top decay yielded by HERWIG and PYTHIA, after 
fitting their hadronization models to LEP and SLD data. Right: Transverse energy of 
the first four jets in W + 4 jet events according to a few generators. 



We wish to present a few studies on the theoretical uncertainty on the top mass 
measurement. Ref. [33] discusses the error in the lepton+jets channel, due to different 
parametrizations of the underlying event in PYTHIA, with the possible inclusion of jet- 
cnergy-scale corrections. As showed in Fig. 4, several tunings and models arc investigated 
and the overall uncertainty is estimated to be Arrit ~ ±0.5 GeV. Scale and parton-density 
uncertainties are instead studied in [34]. Fig. 4 also presents the NLO cross section 
a{pp tt), as a function of nit, for various values of renormalization and factorization 
scales. The relation Amt/mt ~ 0.2 Aa/a + 0.05 emphasizes that the relative error on 
TOt cannot be below 5%, if rrit is determined from a measurement of the cross section. 
Another source of uncertainty on nit is 6-quark fragmentation in top decay, which relies on 
the hadronization models implemented in MC generators, fitted to e"'"e~ data. In Ref. [35] 
it was found that the default parametrizations of cluster and string models are uncapable 
of fitting i?-hadron data from LEP and SLD. After fitting few parameters, PYTHIA 
gives a good description of the data, whereas HERWIG is only marginally consistent. 
Therefore, as shown in Fig. 5, HERWIG and PYTHIA predictions for the B-spcctrum in 
top decay look quite different: such a result will have an impact on the MC uncertainty on 



Table I.: Cross sections in pb for {Z/'y* e+e ) + n jets at the LHC according to 
different codes. The following cuts are set on jets: pt > 20 GeV, \ri\ < 2.5, AR = 
V(A(/))2 + (A77)2 > 0.4. 
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mt- The analysis [36] investigates instead the systematic error due to initial- (ISR), final- 
state radiation (FSR), underlying event (UE) and hadronization on the reconstruction 
of rrit as the invariant mass of the Wb-jet combination, in the lepton-|-jets channel. The 
6-jet is clustered according to the kr (KtJet) and cone (PxCone) algorithms. Ref. [36] 
found that, relying on the kx algorithm, the reconstructed rrit shows visible dependence 
on FSR, ISR and UE, whereas hadronization effects are very little. When using the 
cone one, FSR and hadronization have a large impact, while ISR and UE effects are 
negligible. Given the different features of the two algorithms, using both of them is 
therefore advisable. 

As for matrix- element generators. Table I [37] quotes the results of a few codes on 
Z + n jets, with n < 6 and Z — > e~^e~ , one of the backgrounds for top production. It is 
remarkable that, as long as a process is implemented, the considered codes agree. More 
exclusive studies on W+ jets were carried out in [38], but even there the used programs 
were found to agree. Fig. 5 shows the transverse energy distribution of the first four 
hardest jets yielded by ALPGEN, ARIADNE, HELAC, MadEvent and SHERPA. 

In summary, we reviewed some of the existing event generators for top signals and 
background and presented a few results for the LHC. Given the large numbers of available 
tools, it is always advisable using at least two codes for comparison. 
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